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Neutron Stars in Rastall Gravity 
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We calculate static and spherically symmetric solutions for the Rastall modification of gravity to 
describe Neutron Stars (NS). The key feature of the Rastall gravity is the non-conservation of the 
energy-momentum tensor proportionally to the space-time curvature. Using realistic equations of 
state for the NS interior we place a conservative bound on the non-GR behaviour of the Rastall 
theory which should be < 1% level. This work presents the more stringent constraints on the 
deviations of GR caused by the Rastall proposal. 
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I. INTRODUCTION 

Einstein’s general relativity (GR) is a robust theory for 
explaining the gravitational phenomena. It operates by 
relating directly the curvature of space-time to its energy 
and momentum content by the so called field equations. 
Importantly, its predictions for solar system experiments 
and binary pulsars systems are in perfect agreement with 
the observational data. On larger (cosmological) scales 
we can keep GR as our standard gravitational theory 
by evoking extra unknown components like dark matter 
(DM) and dark energy (DE). 

Alternative gravity theories appear usually when one 
allows for extra degrees of freedom or desires to abandon 
the concept of DM or DE components. Many alternative 
theories have been proposed as scalar-tensor theories (be¬ 
ing Brans-Dicke|lJ a particular case of this class), f(R) 
theories [5] and others J]. Such theories can be tested 
with great precision using solar system tests and binary 
pulsars. For the latter, see [1J. 

In particular, Rastall Gravity can be seen as a gener¬ 
alization of General Relativity. In 1972 P. Rastall 0 
proposed an adjustment to Einstein’s equation which 
causes a violation of the usual conservation law, making 
the covariant divergence of stress-energy tensor propor¬ 
tional the covariant divergence of the curvature scalar, 
i.e., T{? oc R- u . Therefore when the local/global geom¬ 
etry is flat we recover Einstein’s field equations. How¬ 
ever, it seems non-trivial to explain the nature of such 
new source. This can be phenomenologically viewed as 
the appearance of quantum effects in the classical con¬ 
text, as discussed in [6], Interestingly, the issue of non¬ 
conservation of T is a feature also found for instance in 
diffusion models [2 and its connection to modified grav¬ 
ity theories has been discussed in literature in Refs. m- 
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Since the energy-momentum tensor is not conserved 
anymore, there is no variational principle in the context 
of the Riemannian geometry. However, some structures 
similar to that of the Rastall’s theory may be found in 
the context of Weyl geometry [10!. External fields in the 
action may lead also to essentially the same structure as 
in Rastall’s theory (for a discussion of external fields in 
gravity theories, see for example EH). 

Cosmological scenarios based on the Rastall proposal 
(Rastall Cosmology) may be degenerate with the ACDM 
at background and at first order levels, i.e., a viable 
Rastall cosmology can be constructed |T2|. On the 
other hand, astrophysical consequences of Rastall pro¬ 
posal have not yet been widely considered. 

Apart from the cosmological solutions, gravitational 
theories must also provide other solutions for stellar and 
black hole configurations and gravitational wave emis¬ 
sion. Concerning the former, the widely known Tolman- 
Oppeinheimer-Volkoff (TOV) equation allows us to de¬ 
termine the GR static and spherical equilibrium config¬ 
uration of the stellar interior EH El]. 

Stars generate energies by nuclear reactions within the 
stellar interior and gravitational forces are supported via 
radiative convection. This balance can be studied with 
Newtonian mechanics, via the Lane-Emden equation, as 
far as p<Cp (in c = 1 units). Even white dwarfs which 
have the gravitational collapse supported by the electron 
degenerescence, maximal mass ~ 1 .4Mq and a radius of 
order of hundreds of kilometers, are well described by the 
classical theory EH- 

Compact stars have p ss p where the density is approx¬ 
imately the nuclear saturation density. Then the use of 
the relativistic theory, as given by the TOV equations, is 
mandatory. The effects of pressure can not be neglected 
and it is not viable to use of Newtonian theory to describe 
this structure (see jl6| for an alternative viewpoint). 

Since the essence of the Rastall gravity is associated to 
high curvature environments, neutron stars seem to be 
promising laboratories for testing it. Our main intention 
in this work is the establish bounds to deviations to GR 
caused by Rastall’s proposal. 

In the next section we obtain the Rastall version of the 
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TOV equations. Numerical results of such new equations 
for the mass-radius diagram adopting realistic equations 
of state for neutron stars are shown in section [Till Finally 
we conclude in last section. 


II. STATIC AND SPHERICAL SOLUTIONS 

The description of stellar objects is simplified when we 
consider that matter is spherically symmetric distributed 
in a static geometry. The resulting spacetime is described 
by 

ds 2 = B(r)dt 2 — A(r)dr 2 — r 2 {d6 2 + sin 2 9d(j) 2 ). (1) 


The Rastall version of the TOV equation, we call it 
RTOV from now on, can be written in a similar form to 
the original TOV equation as 
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where new effective quantities were defined using the su¬ 
perscript . The new quantities p and p are related to 
p and p by 


The matter content of an ideal fluid, neglecting the mag¬ 
netic field and the rotational frequency, is contained in 
the energy-momentum tensor 

Tfj, v = (p + p)Ufj,u v — pg^v, ( 2 ) 

where u a is the 4-velocity (with = 1), g is the 

component from metric tensor. 

An important step in deriving the standard TOV equa¬ 
tions is the conservation law, 


Tf 1 ; = o, (3) 

which is concomitant with the use of Einstein’s equations, 

R, 1V - ^g^R = 87 tGT^, (4) 


where R is the Ricci tensor and R is the curvature 
scalar (trace of tensor R^). After following these equa¬ 
tions for the geometry defined in 0 we find then the 
TOV equation 

dp GMp (l + ?)(! + 

dr r 2 1 _2GM ’ W 

r 

where on should keep in mind that 

M = f 47t pr' 2 dr\ (6) 

J o 

is a definition used during the derivation of the TOV 
equation and is automatically associated to the mass con¬ 
tained within r. See mm for a revision of this calcu¬ 
lation. 

Let us now redo the same steps in finding ([5]), but using 
however the Rastall gravity. The conservation law will be 
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where the GR framework is recovered by setting A =1. 

One can recast the field equation 0 according to the 
Rastall’s proposal [5] as 
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p = a\p + 3 a 2 p, 
p = a 2 p + a 3 p , 


where 
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Here we assumed that A = 1 + 77 and 77 is the new term 
which parameterizes deviations from General Relativity. 

The signs of 07 , 0:2 and 03 are related to the stellar 
stability by equation 0. Firstly, note that the effective 
Rastall mass M depends now on p and p. Also, if p < 0 
the equilibrium in equation 0 is lost. Therefore, from 
( 11 ), if the quantity 3 p dominates the density p then 


when 02 < 0 leads to p < 0 and M < 0. Otherwise 
if the density p dominates then p < 0 and there is no 
equilibrium. Another aspect which we should note from 
the definitions ( 12 ) is that 01,02 and 03 are negative in 


the following interval 77 = (—2/3,—1/2), 77 = (—1/2,0), 
and 77 = (—2,—1/2), respectively. Consequently, values 
77 < 0 may be problematic since this leads to 02 < 0 . 
Remark that p + p = p + p and the null (and also the 
strong) energy condition seem to be preserved, but under 
extreme conditions 02 negative may lead to p < 0 , im¬ 
plying violation of the energy conditions and generating 
instabilities. The numerical calculations to be presented 
in the next section confirm this possibility. 


III. NUMERICAL SOLUTION 

For solving the equilibrium equation (|9j) we must deter¬ 
mine the EoS of the stellar interior which, in general, has 
as a free parameter the central density (p(r = 0) = po). 
For a giving central density value we obtain the corre¬ 
sponding mass and a radius of the star, i.e., one single 
point in the mass-radius diagram. One should vary (po) 
in order to obtain a curve in this plane. This is the so 
called mass-radius diagram. 

The inner structure of the NS is still not well under¬ 
stood. Many groups have developed possible EoSs us¬ 
ing microscopic many-body calculations based on, for in¬ 
stance, phenomenological relativistic mean-field theory 
and nucleon-nucleon interactions mm- 
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In order to get an insight on the internal structure of 
neutron stars let us firstly assume a simple pure neu¬ 
tron star with nucleon-nucleon interaction. In this case, 
a polytrope can be used as an approximation for the neu¬ 
tron matter equation of state 

p = up 1 . (13) 

Following the description of an interacting Fermi gas 
model EoS as in Ref. [22] we have n ~ 2.0 x 10 5 
when 7 = 2 . 

The corresponding mass-radius diagram is shown in 
figure 0 Since the TOV equation is a non-trivial 
integro-differential equation we have used an adaptive 
step-size Runge-Kutta method as described in [23]. 

Note that we recall in this plot the usual definition for 
mass as presented in 0 . The labels in the figure indi¬ 
cate the 77 values. The black line represents the standard 
GR case ( 77 = 0). For 77 = 10 -3 the mass-radius pro¬ 
file around the maximum mass is very close to the GR. 
Large radii equilibrium configurations in the Rastall case 
are more massive than in GR. The red line (77 = 0.05) 
clearly fails in predicting the existence of typical neutron 
stars with masses around 2M 0 within a 10 — 15 km ra¬ 
dius. It is worth noting that the numerical evaluation 
does not converge for 77 < 0. In this case the hydrostatic 
equilibrium equations |9| are unstable and therefore no 
stable stellar configurations are found. This indicates a 
possible violation of the energy conditions. 



FIG. 1: Mass-radius contours for a pure neutron star (with 
interactions) using Prakash Method. 

Let us concentrate now our discussion on more realis¬ 
tic EoSs. The BSk family of equations aim to provide 
an unified description of the stellar interior treating con¬ 
sistently transitions between outer and inner crust (and 


core). They have 23 free coefficients which have to be fit¬ 
ted numerically. Each BSk equation of state correspond 
to one specific numerical fit. At neutron star densities, 
the unified BSkl9, BSk20, and BSk21 EoSs approximate, 
respectively, the soft, moderate, and stiff astrophysical 
EoSs FPS [Ml |25!, APR [26], and V18 [23 (see also [28] 
for discussion and references). 

Solving now the differential equation ([ 9 ]) for these EoSs 
and then translating our results to the usual definition 
of mass 0 as before, we find the mass-radius diagrams 
shown in Figs. 0, 0 and Q, respectively. 

Particularly, the BSk models have slightly distinct pre¬ 
dictions for the maximum mass. In general, BSkl9 EoS 
does not allow masses larger than 2M 0 . Then, observa¬ 
tions of very massive neutron stars with a radius ~ 13 km 
would favor the BSk20 and BSk21 fits. Masses values of 
some observed neutron stars are usually compiled in [23 ■ 

Both EoSs place a similar bound of the Rastall param¬ 
eter. Clearly, 77 = 0.05 leads to unrealistic configurations. 
Positive 77 values showed here slightly change the maxi¬ 
mum mass of the star. However, the most sensitive effect 
of changing 77 values is to shift equilibrium configurations 
to larger radii. For example, 1.5M 0 objects tipically have 
radii around 10 — 12 km in GR (77 = 0). If 77 = 0.01 these 
radii values are shifted to the range 14 — 17 km. Cur¬ 
rent uncertainties on the radius determination allow us 
to barely differentiate neutron stars radius at this level. 
However, it is quite useful to quote here recent results 
on the determination of R presented in the literature. 
In Ref. [30], for example, bayesian statistics has been 
employed to deal with the free parameters involved in 
modeling neutron stars. It has been found that, for the 
preferred equation of state, the radius of a 1.5 M 0 object 
is R 1.5 = IO.SIq ; 3 & 171 and in general the peak probabil¬ 
ities are highly clustered in the 9 — 12 km range. Also, 
taking into account the results reviewed in Ref. m we 
can state that neutron star radii in the range 14 — 17 
km should be severely disfavored. On the other hand, 
there are also different estimates of R which are not so 
restricted leaving some room for admitting neutron stars 
with radius around R 1.5 ~ 15 km [521153] - Importantly, 
the only safe claim concerning the observed radius of neu¬ 
tron stars is that R 1.5 > 15 km are severely disfavored. 
The difference among the predicted values of radii occurs 
particularly because of the largely unknown systematic 
bias introduced by using the blackbody model instead of 
an atmosphere model or by a priori assumptions on the 
chemical composition of the atmosphere (see [H] for a 
review). 

From the above considerations It is worth noting that 
giving the current uncertainties in the determination of 
neutron star masses and radius, values of order 77 < 0.01 
provide stellar configurations which are consistent ob¬ 
servations and almost indistinguishable from GR predic¬ 
tions. 
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FIG. 2: Mass-Radius contours for representative values of 77 
in the BSkl9 model. 
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FIG. 3: Mass-Radius contours for representative values of 77 
in the BSk20 model. 

IV. CONCLUSION 

In this text, we have discussed the fate of stellar equi¬ 
librium in the context of the Rastall’s theory of grav¬ 
ity. Essentially, we have considered the case of neutron 
stars, where the matter conditions are extreme and the 
deviations of General Relativity may lead to some im- 
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FIG. 4: Mass-Radius contours for representative values of 77 
in the BSk21 model. 

portant effects. In fact, using three equations of state 
for the neutron stars, from the family of BSk models, we 
find that only small departures of the General Relativity 
are consistent. These strong restrictions on the violation 
of the energy-momentum tensor conservation laws may 
break the degeneracy that appears at cosmological level 
between Rastall’s theory and the standard cosmological 
model based on General Relativity, the ACDM model 

m- 

Since strong field effects are present in stellar systems 
like neutron stars, they seem to be an excellent labo¬ 
ratory to test deviations from the standard theories of 
gravity, in spite of the uncertainties existing in the study 
of the nuclear matter configuration in extreme density 
regimes. The obtained restrictions on the Rastall’s pa¬ 
rameter (of the order of < 1% with respect to the General 
Relativity value, A = 1) are, on the other hand, consistent 
with the interpretation of such deviations as manifesta¬ 
tion of quantum effects in gravitational systems. Precise 
future observations of R 1.5 can improve our bound to the 
< 0 . 1 % level. 

Numerical studies, and considerations based on the en¬ 
ergy conditions, restrict the Rastall’s parameter to A > 1. 
Curiously, in the study of cosmic strings in the frame¬ 
work of Rastall’s theory carried out in reference [35] it 
was found that the space-time of the cosmic string with 
abelian fields becomes singular for A > 1, i.e., 77 > 0. 
However, this apparent contradiction must be faced with 
care. The absence of a Lagrangian in Rastall’s gravity 
makes the full determination of the equations of mo¬ 
tion for the matter fields challenging. In the reference 
[35] , only the minimal extension of the General Relativ- 
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ity case was taken into account. To some criticisms on 
the Rastall’s proposal based on hydrodynamics, see ref¬ 
erence 1551 . 

Finally, concerning the bounds we have proposed for 
77 in the Rastall’s theory, there is a remark to be made. 
The Rastall’s theory leads, in the weak field limit, to a 
Poisson equation but for an effective gravitational cou¬ 
pling that is related to the usual gravitational constant 
G by terms of order rj. Measurements of the gravita¬ 
tional constant have a precision of order of 10 - 4 m, 
but use local experiments like torsion balances or atomic 
interferometry |38| . Even if our results may somehow be 
degenerated with the uncertainties on the knowledge of 
the gravitational coupling the uncertainty in measuring 
G can not affect the present result at the level of 1%. 
There is also a degeneracy with the uncertainties coming 
from the knowledge of the equation of state of neutron 
stars. Then more reliable equations of state and, espe¬ 


cially, more precise measurements of neutron star radii 
are crucial to break this degeneracy. In this sense, in 
the present state of art, the results obtained may indi¬ 
cate that Rastall’s theory is still compatible with obser¬ 
vations (in what concerns neutron stars) for a small range 
of values of A, but this range can become even narrower 
in the future. The results obtained in this present work 
reinforce that the high-density environments structures, 
like stars, may constraint very strongly deviations from 
General Relativity and, in particular, deviations from a 
possible violation of the usual conservation laws. 
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